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ABSTRACT 

We investigate seed magnetic field generation in the early universe by radiation force of first stars. In 
a previous study with the steady assumption, large amplitudes(~ 10 -15 G for first stars, ~ 10~ n G for 
QSOs) are predicted. In this study, we formulate this issue in an unsteady framework. Then, we consider 
a specific model of magnetic field generation around a very massive first star. Consequently, we (1) find 
that the steady assumption is not valid in realistic situations, and (2) obtain much smaller magnetic 
field strength than predicted by Langer et al. In addition, we find that the momentum transfer process 
during photoionization is more important than Thomson scattering. The resultant magnetic flux density 
around the first star is < 1CP 19 G. This seed magnetic field will not affect subsequent star formation in 
the neighborhood of first stars. 

Subject headings: early universe — HH regions — magnetic fields — radiative transfer 
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1. INTRODUCTION 



One of the most important current issues in cosmology is understanding the star formation history of the universe, 
especially the early universe. Owing to progress in the last decade, mainly due to theoretical studies, we now expect 
the very first stars in the early universe to be very massive ( > IOOMq), compared to those in present-d ay galaxies(e.g. 

|2002 |Nakamura fc Umemura|2001| Abel, Bryan, Norman] 



Bromm, Coppi & Larson 



2002 Yoshida 



2006) 



The chief reason that the first stars are very massive is the lack of heavy elements in the primordial gas, since the heavy 
elements are efficient coolants. Theoretical studies suggest that when the metallicity of the star forming gas exceeds 
10 _5 tolO _6 Zf7i, the star form ation mode transits from very massive star formation to normal star formation^ Schneider et 



-0: 



al.|2003 Omukai et al.|2005 |. Thus, metallicity is the key factor for the transition of the star formation mode7 

nother possible parameter is magnetic field strength. In the present-day star-forming molecular cloud, energy density 
of the magnetic field is comparable to gravitational or turbulent energy density, while it is assumed that the magnetic 
field is weak in the primordial gas. Magnetic field in present-day galaxies basically suppresses the star formation process. 
The effects of magnetic field on the early star formation episode in the universe hav e been studied by several authors. 
The coupling between magnetic field and pure primordial gas has been investigated( |Maki fe Susa 2004 2007). They 
found that the magnetic field is always frozen in the gas in so far as the magnetic field strength is not strong enough to 
prevent gravitational contraction. Based on the frozen-in assumption, it is found that bipolar outflows emerge in cases 
B > 10 _9 G is satisfied at nn = 10 3 cm~ 3 ( Machida et al.|[2006 ). Tan & Blackman (2004) suggest that magnetorotational 
Instability (MRI) is activated in the accretion disk surrounding proto first stars in cases where B > 10~ 10 G at tih 



10 3 cm 3 , which results in efficient angular momentum transfer in the disk. It is also reported by Schleicher et al. 



(2009) 



that gravitational contraction of the primordial gas is delayed in cases where the comoving field strength is larger than 
10- 10 G. 

Therefore, it is crucial to determine the magnetic field strength that exists in the environment of the star-forming gas 
cloud in the early universe in order to understand the transition of the star formation mode. Magnetic field strength 
at the star for mation site in the early universe is basically provided by the field generated prior to no nlinear structure 



formation (e.g., Harrison 1970 Mishustin fc Ruzmaikin p.972; Turnerfc Widrow 1988 Ichiki et al 



field generat ed by t he Bicrmann battery effect (Biermann 1950) during the structure formation (e.g. 



2006), as w ell as the 
Xu et al 



2008|) 



or 

°G, 



reionization ( Gnedin et al. 2000 ) . The predicted field strength in these models is within the range of 10 
which is too weak to affect star formation. Another possibility is that the magnetic field spewed out from the formed 

A G like normal stars, a simple 
10~ 16 G at intergalactic medium(IGM) densities. 



stars (e.g., Bisnovatyi-Kogan et al.|1973 [). In cases where the field strength is as large as 
flux freezing assumption again leads to 



On the other hand, |Langer et al. ( 2003 ) suggest a mechanism to generate a large amplitude of seed magnetic field as 



large as 10 _11 G, which could affect the star formation in the early universe. They consider HII regions formed by luminous 
sources like QSOs or first stars. The radiation force field due to the photons emitted by these objects is anisotropic because 
of the dumpiness of surrounding media. This anisotropy causes charge separation in the HII region, which generates a 
large scale eddy current as well as the magnetic fields. However, the assumption in this theory is idealized. In their 
theory, anisotropy of radiation is dealt with as perturbation, and the HH region is assumed to be steady. In reality, dense 
shadows are generated behind the nearby dense halos. In addition, HE regions formed by QSOs or first stars are actually 
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unsteady, a nd ionization fronts never reach the Stromgren radius ( Shapiro fc Giroux|1987 Kitayama et al.|2004 Yoshida 
eFll][2007| 



In this paper, we investigate t his magnetic field gen eration process by radiation force in nonlinear and unsteady frame- 



work along the lines of theory of Langer et al. (2003). We consider a specific model of magnetic field generation around 
a very massive first star and assess whether the generated magnetic field affects the subsequent star formation. 

In Section 2, we describe the basic equations. A brief description and the setup of our models are described in Section 
3. In Section 4, we show the results of our calculation, and Section 5 is devoted to discussion. Finally, in Section 6, we 
conclude our research. 

2. BASIC EQUATIONS 

2.1. Magnetic field generation 

In this section, we describe the basic equations of magnetic field generation. To begin with, the force balance on a 
single electron is written as follows: 

Q = -eE+-j + f l . &d - V ^--v c xB (1) 

c c n c c 

Here, e, E, B, and j denote the elementary charge, electric field strength, magnetic flux density, and electric current 
density, respectively. p B , n c , and v c represent the pressure, number density and velocity of electrons as fluid, respectively. 
<7 C is the conductivity and / rad is the radiation force on an electron. 
The fluid velocity v is defined as 

ripiripVp + n e m e v e . . 

v = p ~ Up 2 

n p m p + n m e 

where n p ,m p ,m c , and v p denote the number density of protons, the mass of protons, the mass of electrons, and the 
velocity of protons, respectively. Thus, the velocity of electrons is approximated as v c ~ v — j/(en e ). Substituting this 
notation into EquationpJ) , we obtain 

Q = -eE+ -j + / rad - — ^ - -v x B + jxB (3) 

<r c n c c cn e 

Since we consider the expanding HII region around a first star, we have to solve the following photoionization rate 
equation for electrons. 

dn t 

~m 

Combining Equation Q with the similar equation for protons, we obtain the charge conservation equation 

| + V.,- = 0, (5) 

where charge density p and current density j are related to the number density of electrons and protons as p = e (n p — n c ) 
and j = e (n p v p — n c v c ). We also need the set of Maxwell equations: 

V-E = 4np, Vx £ = -- — , V • B = 0, VxB=^j + ^- (6) 

Equations (|3|,Q,((5| and (pi) combined with hydrodynamic equations are the fundamental equations for this problem. 

In order to obtain the equation that describes the growth of magnetic flux density, we apply — (47rcr c /ec)Vx to both 
sides of Equation ^. We have 

47r . 47rer r „ 47rav . „. 47rtr r ( j x B\ Aiio~c Ana c 

= V x j + ^V xE+ -^V x(vxB ^V x J - |Vn c x V Pc ^V x / rad 

c c c z c z \ n c J cen£ ce 

Using the second and fourth equations of Maxwell Equations (p| , we have 

1 d 2 B -jz-d Airo- c dB Aira c 47rcr c f j x B\ 4ttct c 4tt<7, 



V • (n e v e ) = T - an e n p (4) 



, -V'B + —^— = — ^V x(vxB)- ^V x J - V« c x V Pc ^V x / rad (7) 

c z at* c z at c z c z \ en c J cen£ ce 

In this expression, Aira c denotes the inverse of the timescale that describes t he converge nce of electron velocity to the 
terminal velocity. This timescale is quite short, since a c ~ 6.5 x 10 6 r 3 / 2 e.s.u.( Lang||l999| ). It is as short as a picosecond 
in cases where we consider photoionized gas. Thus, the first and second terms on the left-hand side of the above equation 
are negligible compared to the third term, since the typical timescale is the lifetime of the source star which is as long as 
Myr, and the typical length scale is kpc. Finally, we obtain 

dB c c 

-r-=Vx(vxB) -Vn e xVft--Vx/ rad (8) 

at en£ e 

Here we also omitted the Hall current term which is proportional to j x B, since this is a higher order term than 
the others, in cases where we consider the generation of very weak magnetic field. The first term on the right-hand side 



denotes the advection term. If the second and third t erms are neglect ed, the frozen-in condition is satisfied. The second 
term describes the well-known Biermann battery term(Biermann|1950 ), and the third term represents the radiation term. 



Thus, the non-zero rotation of the radiation force field always introduces the source term of magnetic field generation. It 
is worth noting that the gravitational force cannot be the source of magnetic field generation, since it is described by the 
gradient of a scalar potential. 

We also mention that the final Equation (pi does not contain the terms directly originated from the electric current. In 
fact, Equation ([8]) is independent of conductivity, o~ c . The time derivative of magnetic field simply comes from Faraday's 
induction equation, i.e., the second equation of Maxwell Equation (km). Thus, the magnetic field generation by the radiation 
force in this paper can be regarded as the consequence of electromagnetic induction effects. The electric field with non-zero 
rotation is generated so as to cancel the radiation force. Thus, it is very natural that the second equation of Maxwell 
Equation (p| coincides with Equation ([8]), if we assume eE = / rad and neglect the advection term and Biermann term 
in Equation (|8| . We also note that a similar problem was considered in deta il for the problem of fiel d generation due to 
radiation forces in accretion disks around black holes (Bisnovatyi-Kogan & Blinniko v|1977|[Contopoulos fc Kazanas|1998| 
Bisnovatyi-Kogan et al.|[2002 Contopoulos et al.|2006||Christodoulou_ et al.|2008p . 



Aside from the magnetic field generation described by Equation d8j), we can assess the charge separation which drives 
the electric current. Taking the divergence of Equation pi) combinea with the charge conservation Equation (pi) and the 
first Maxwell Equation®, we have 

% + 4wa c p = a c V-(--vxB + —j xB-^ + -/rad ) (9) 

ot \ c en c c en c e / 

The second term on the left hand-side is much larger than the first term, since l/(4%a c ) <C t is always satisfied. Therefore, 
the charge density is assessed as 

P=]-V-(--vxB+— jxB-^ + -f iad ) (10) 

4ir \ c en c c en c e / 

Thus, separated charge is piled up where the divergence of the force field is large. For instance, if we consider only / rad 



on the right-hand side of Equation ( 10 ), V • / rad becomes negative in the neighborhood of the ionization front, since / rad 



changes very rapidly across the ionization front. Thus, electrons will accumulate around the ionization front. 

2.2. Momentum transfer from radiation to electrons 

In order to follow the generation of magnetic field, we need to evaluate the radiation force on electrons, / rad . Two 
elementary processes could be potentially important for this momentum transfer. First of all, Thomson scattering could 
contribute to / rad . We use the formal solution of the radiation transfer equation for / rad T as follows: 

er T f UL a T f°° 

/rad t = — / F 0v di/-\ / F 0l/ eaq?[-T VL a(v)]dv, (11) 

c Jo c J VL 

where / rad T denotes the radiation force per single electron due to Thomson scattering, <tt is the cross section of Thomson 
scattering, Fq v is the unabsorbed energy flux density, vj, denotes the Lyman-limit frequency, r Vh is the optical depth at 
Lyman limit regarding the photoionization, and a(y) denotes the frequency dependence of photoionization cross section, 
which is normalized at Lyman limit, i.e., a(fi,) = 1 is satisfied. 

Second, the photoionization process itself also transfers the photon momentum to electrons. Its contribution to / rad is 
given as 

/rad.i=T— / <Tv L a(y)F 0v exp[-T VL a(v)]dv, (12) 

2 cn c J VL 

where a VL denote s th e photoionization cross section at the Lyman limit. We remark that the factor 1/2 in the right-hand 



side of Equation ( 12 ) is due to the fact that the photon momentum is equally delivered to protons and electrons. In fact, 
when electrons and protons are emitted by photoionization, its momentum is basically delivered to surrounding electrons 
and protons by Coulomb interaction. Since the number density and charge of each particle is the same for electrons and 
protons, electrons gain half of the momentum transferred from photons. Because of the large inertia of protons, we can 
also safely assume that only electrons are accelerated by this momentum transfer process. 

The photoionization rate T in Equation Q is also obtained by assuming the formal solution of the radiation transfer 
equation: 

<J VL a(v)-^-exjp[-T VL a{v)]dv (13) 



In practice, the photon conserving scheme (Abel, Norman fc Madau|1999 Susa|2006) is used when we solve Equation Q 



so that we can trace the propagation of the ionization front properly. 

3. THE MODEL 

3.1. Description of the model 
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Fig. 1. — Schematic view of the magnetic field generation process. Dashed arrows represent the direction of electron flow. 

We consider gas surrounding a first star. This ambient gas is ionized by the radiation from the source, forming an HII 
region. The momentum transfer from ionizing photons to electrons leads to a slight charge separation in the HII region. 
It creates a weak electric field balancing the radiation force. 

The electric field and the distribution of charge cannot be isotropic, because the propagation speed of the ionization 
front is anisotropic due to the dumpiness of the surrounding media. As illustrated in FigureJT] overdense regions in the 
neighborhood of the radiation source form shadows behind themselves that result in anisotropic radiation flux as well 
as the electric field/charge distribution. Consequently, a transverse electric field along the ionization front is generated. 
This field could drive eddy currents surrounding the boundary between the shadowed and unshadowed regions. Thus, the 
magnetic field is expected to be generated at this boundary region. 

3.2. Setup 
We consider a first star of 500Af©(L = 5 x 10 40 ergs" 1 , t agc = 2 x 10 6 yr, T g = 10 5 K (|Schaerer | |2002[) ), at redshift 



z ~ 20. Since the ionization front breaks out the small host halos with Mhaio ^ 10 7 MQ( Kitayama et al.|2(J(J4[ ) , we consider 
the expansion of the HII region in intergalactic space. The initial number density of the ambient gas in intergalactic space 
is no = 10~ 2 cm~ 3 at such redshift. It is also assumed that the gas is neutral when the source star is turned on. 

As shown in Figurd2l we follow the growth of magnetic flux density in a two-dimensional model. The computational 
domain is a strip of 500pcx5kpc. The radiation flux flows from the left along the long edge of the domain. The assumed 
flux on the left edge is the flux lOpc away from the source first star. We set an overdense region of nocm -3 , Z?kpc 
away from the left edge of the simulated region, which casts a shadow behind it (Figurc[2]). We tested four models 
A:n Q = lcm- 3 /D = 2kpc, B:n = lcm" 3 /!) = 200pc, C:n = lOcm" 3 /- = 200pc, and D:n = lOcm- 3 /^ = 2k P c - 
The overdense region has a core-halo structure. The diameter of the core (n = nocm -3 ) is 50pc, which is surrounded by 
the envelope of n oc r~ 2 , where r is the radial distance from the center of the overdense region. 

In this paper, we assume that the gas is isothermal (T — 10 4 K), which results in eliminating the Biermann battery 
effect. In addition, we also assume that the gas is at rest, i.e. v = 0. This assumption is valid when the ionization 
front is propagating in the intergalactic space; however, it is not appropriate in cases where the dense gas clump is 
photoevaporated because of the thermal pressure of the photoheated gas. These effects are left for future work. In this 
paper, we focus on the effects of radiation. 



We also remark that the timescale for the ionization front to reach the Stromgren radius is l/(cm e ) ~ 1.2 x 10 7 yr(Shu 



1992), which is sufficiently longer than the lifetime of the source star. T herefore, the HII region does not settle down to 



a steady state during the lifetime of the source star( |Yoshida et al.||2007[ ). 

4. RESULTS 

4.1. Orders of magnitude 

Before we move to the numerical results, it is useful to show the expected magnetic field strength by rough calculations. 
Following the isothermal and static fluid assumptions, Equation ([8]) is reduced to a very simple equation: 

f =-;vx/, d4) 

This expression shows that the magnetic field is effectively generated where the radiation force field has a large shear. 
It is obvious that / ra( j should have a large derivative at the boundary between the shadowed and unshadowed regions. 
Thus, we focus our attention on such boundary region. 



The radiation force / rad is assessed by Equations (11) and (12) as follows: 

"HI 



/rad — /rad.T + / 



rad. I 



' V 2n c 

where F is the total energy flux coming from the source star. Thus, B is evaluated as 
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x 10 ergs 
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(15) 



(16) 



Viopcy V 5 ' "°" J \zkpcj V 2xl ° 6 y r / 

where Ar denotes the length scale that represents the sharpness of the boundary between the shadowed and unshadowed 
regions. Here we assume nm/n ~ 1 and the flux is not absorbed significantly (F ~ L*/(47ri? 2 )). This assumption is only 
valid in the neighborhood of the ionization front. I n t he neutral region, dumping by the absorption of ionizing radiation 
reduces the radiation force. Therefore, Equation ( 16 ) should be regarded as an upper limit of the generated magnetic 
field strength. 

It is also worth noting that Thomson scattering is negligible in neutral regions since the cross section is smaller than 
photoionization by 7 orders of magnitude, although it could be dominant in the highly ionized region. That is why the 
contribution from / radT is omitted in Equation (16). 



4.2. Numerical results 



We integrate Equation (14) through the lifetime of the source star (2 x 10 6 yr) after it turned on. Three snapshots 
for the distribution of magnetic field strength are shown in FigurtkMor model A:n c = 1cm -3 , D = 2kpc. Three panels 
correspond to 0.5, 1 and 2Myr after the ignition of the source star, from top to bottom , respectively. The radiation 
propagates from left to right in the panels, as the ionization front(green curve). The ionization front is trapped at the 
overdense core located 2kpc away from the left edge. Note that the right panels are compressed along the direction of the 
light rays, since the shadow is a quite elongated structure. The left panels have the same scale in both the horizontal and 
vertical directions. 

The orange area denotes the region where the magnetic field is perpendicular to the plane of the page and directed to its 
front side. In the blue region, the magnetic field points in the opposite direction. First of all, it is clear that the magnetic 
field strength peaks at the boundary of the shadow cast by the dense core. As discussed in Section |4~Tj / rad should vary 
significantly across the boundary between the shadowed and unshadowed regions, since the radiation flux is very weak in 
the shadow. In fact, as shown in FigureHJ where the color contour of / ra( j is shown, / rad changes significantly behind the 
dense core. Since the ionization front is trapped at the dense core, the large shear of / rad is kept relatively long, while it 
is only transient in the low-density region far from the dense core. Thus, the magnetic flux density reaches its peak just 
behind the dense core, since it has enough time to grow. 

We also remark that the border of the shadow does not coincide with the ionization front. As shown in Figurd4l / rad 
is larger behind the ionization front, i.e., in the neutral region. This behavior indicates that the momentum transfer from 
photons to electrons through the photoionization process is more important than Thomson scattering since / rad j is larger 
in the neutral region, while Thomson scattering is important in the ionized region. 

The time evolution of the peak magnetic field strength is plotted in Figur^EJ The field strength immediately reaches 
the level of ~ 10" 21 G after the ionization front hits the overdense region at t = 2 x 10 5 yr. After that, the field strength 
grows gradually until the source star dies. 

We also show three snapshots at 2Myr for three other models in Figure[6] These panels correspond to models B, C and 
D from top to bottom, respectively. 
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Fig. 2. — Schematic view of the computational domain. The shadowed region is formed behind the overdense region , while the other region 
is exposed to the radiation field. 
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Fig. 3. — Color contour maps of the magnetic field strength at three snapshots at 0.5, 1 and 2Myr. The three snapshots on the left are 
extended views around the dense core. Green lines denote the position of the ionization front. Orange represents the region in which B is 
directed to the front side of the page, whereas blue is used for the opposite direction. 
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Fig. 4. — Color contour map of the radiation force / ra( j lMyr after the ignition of the source star. White arrows schematically show the 
flow of electrons. 
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Fig. 5. — Time evolution of the peak magnetic field strength is plotted for model A:n c = 1cm , D = 2kpc. 
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Fig. 6. — Color contour maps of the magnetic field strength at 2Myr for model B:n c = 1cm -3 , D = 200pc; model C:n c = 10cm -3 , 
D = 2kpc; and model D:n c = 10cm -3 , D = 200pc. The panels on the left represent extended views around the dense core. Green lines 
denote the position of the ionization front. Orange represents the region in which B is directed to the front side of the page, whereas blue is 
used for the opposite direction. 



In model B the dense core is located 10 times closer to the left edge than in the case of model A. We expect 100 times 
stronger magnetic field strength following the scaling relation of Equation ( 16 ). However, we have the maximum magnetic 
field strength of ~ 2.0 x 10 _20 G. In this case, the ionization front is not trapped by the dense core. Consequently, the 
magnetic field does not have enough time to grow since the shadowed region behind the dense core is feathered after the 
ionization front passes through. Thus, the field strength does not follow the simple scaling relation. 

Model C corresponds to the case where the density of the core is 10 times larger than that in model A, while the 
distance D is the same. The peak magnetic field strength in this model is as large as 3.0 x 10~ 21 G, which is smaller 
than the case of model A by a factor of few. Since the ionization front does not pass through even at the envelope of the 
overdense region, the rotation of / rad (i.e., the shear of / rad ) around the dense core becomes relatively smaller than that 
in model A. Therefore, the magnetic field strength is smaller than the case of model A. 

Finally, we describe the results of model D, which are shown in the bottom panels of FigurelS] In model D, the core is 
located 10 times closer, and is 10 times denser than that in model A. As shown in FigurelHJ the ionization front is trapped 
in front of the dense core because of the high enough core density. It is also clear that the ionizing radiation flux is high 
enough to ionize the envelope, unlike model C. In this case, the maximal magnetic field strength is ~ 10 _19 G. 

The parameters employed in models A to D scan the reasonable range. In fact, the actual size of the HII region of a 
500M Q first star is as large as several kpc, and the size of the first halos is < lOOpc. The density of the virialized halo 
is ~ 1cm -3 at z = 20. Thus, we can conclude that the magnetic flux density generated by radiation from first stars is 
<10- 19 G. 



5. DISCUSSION 

We found that the generated magnetic field strength is ~ 10 _19 G at maximum if we ad opt reasonable ran ges for the 
parameters in our model. This value is much smaller than those predicted by steady theory) Langer et al.|2003 ), although 
we take into account the additional contribution by / rad j. If we consider only Thomson scattering as done in Langer et 
al. (2003), we obtain a magnetic field of ~ 10 _25 G. This discrepancy comes from the steady assumption. In their study, 
the magnetic field strength was estimated using steady equations, however; as shown in Section |2.1[ a much longer time 
than the age of the universe is needed to settle down to steady state in this situation. Considering a realistic timescale, 
such as the age of the source star ~ 10 6 yr, we conclude that the steady assumption is wrong. 

In the present work, the system is assumed to be isothermal (= 10 4 K) for simplicity, which results in eliminating the 
Biermann battery term. In a more realistic situation, we have to consider the effects of radiation hydrodyna mics as well 
as th e heating and cooling of gas by which we can compute the magnetic field generated by the battery effect) Biermann 
1950). In fact, the order of magnitude of magnetic field generated by the battery effect is 
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/ Ar 
\10pc, 

where Ar denotes the typical length of n c and p c change significantly and 9 denotes the angle between Vn c and VT. 
This expected field strength is even larger than the flux density due to radiation pressure, although we have significant 
uncertainties in the above estimate. Thus, we need to include the Biermann battery term to assess the correct magnetic 
field strength utilizing radiation hydrodynamic calculations, which is left for future work. 

We also remark that the magnetic field strength obtained in this paper seems insufficient to affect star formation in the 
early universe, since B ~ 10 _19 G at IGM cannot account for 10~ 10 tol0~ 9 G at n = 10 3 cm -3 required for jet formation 
(Machida et al. 2006|) and MRI activation) Tan & Blackman 2004). However, it is also worth n oting that a very re cent 



study by Schleicher et al.| (2010) indicates that turbulence in the host halos of first stars (e.g., Wise & Abel 2008) can 
amplify the magnetic field to equipartition strength ve ry quickly. If this is t rue, this might suppress the fragmentation of 
gas disks formed at the center of the star- forming gas ( Machida et al.| [2008), and help the MRI activation. 



6. CONCLUSION 



We investigate the ma gnetic field generatio n process by radiation force in a nonlinear and unsteady framework along 



the line of the theory of Langer et al. ( 2003 ) . We consider a specific model of magnetic field generation around a very 
massive first star. Using analytical and numerical calculations, we find that the steady assumption is not valid in realistic 
situations and obtain a much smaller magnetic field strength than predicted by the steady theory. In addition, we find 
that the momentum transfer process during photoionization is more important than Thomson scattering. The resultant 
magnetic flux density around first stars is ~ 10 _19 G. This seed magnetic field does not seem to affect subsequent star 
formation in the neighborhood of first stars. 
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